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A summary  is  given  of  leehniques  and  results  for  calculating  light  scatter- 
ing in  anisotropic  crystals  valid  lor  arbitrary  directions  of  the  incident  beam, 
the  scattered  beam  and  the  crystal  surface  normal  relative  to  the  crystal  axes. 
A dyadic  Green's  function  that  distinguishes  ray  and  propagation  directions 
leads  to  a scattering  efficiency  inside  the  crystal  that  involves  the  Gaussian  cur- 
vature of  the  surface  of  wave  normals.  After  taking  account  of  the  solid  angle 
expansion,  and  the  source  volume  demagnificalion  at  the  crystal  surface,  a 
scattering  formula  is  given  suitable  for  comparison  with  experiments  done  out- 
side the  crystal.  Application  is  made  to  Brillouin  and  Raman  scattering. 


I.  INTRODUCTION 

The  availability  of  laser  sources  has  led  to  a great  upsurge  in  light  scattering  experiments 
in  recent  years.  Many  experiments  have  been  performed  on  crystals  of  low  symmetry  and 
large  optical  anisotropy  while  a few  experiments 12  have  used  light  scattering  to  measure  quan- 
titatively the  nonlinear  polarization  causing  the  scattering.  Such  numerical  measurements  re- 
quire a theory  which  (l)  accounts  correctly  for  the  noncollinenrity  of  the  ray  vector  and  the 
wavevector  that  results  from  the  optical  anisotropy  and  (2)  relates  the  nonlinear  polarization  to 
measuring  instrument  parameters  which  are  outside  the  crystal.  Previous  theories3  have  often 
mishandled  the  first  requirement  and  have  always  ignored  the  second. 

The  purpose  of  our  study  of  the  theory  of  light  scattering,  of  which  we  present  a sum- 
mary here,  has  been  to  remedy  these  deficiencies  of  previous  theories.  To  remedy  the  first 
deficiency  resulting  from  optical  anisotropy  we  have  based  our  theory  on  a Green’s  function 
solution  of  the  electric  field  wave  equation. 4,5  This  technique  is  appropriate  because  the 
scattering  volume  is  typically  very  small  compared  to  the  size  of  the  crystal  under  study  when 
laser  sources  are  used.  Because  even  the  crystal  surface  is  in  the  far  field,  an  asymptotic 
evaluation  of  the  Green's  function  inside  the  crystal  is  adequate.  Our  first  - and  incorrect  - 
evaluation4  used  a stationary  phase  integration  over  two  variables  followed  by  a residue  in- 
tegration over  the  third  variable.  The  appealing  but  nevertheless  incorrect  result  was  what 
would  have  been  calculated  for  an  isotropic  medium  with  the  subsequent  replacement  of  the 
isotropic  refractive  index  by  the  one  appropriate  to  the  correct  ray  direction  in  the  crystal. 
This  replacement  procedure  has  been  used  by  others.6  Our  second  - and  correct  - evaluation5,7 
follows  that  of  Kogelniks  and  Kogelnik  and  Motz'1  who  used  it  on  a magnetoionic  media  prob- 
lem. Their  work  was  based  on  a stationary  phase  technique  of  Lighthill10  in  which  the  residue 
integration  was  performed  first  and  the  stationary  phase  method  was  applied  second.  The 
correct  asymptotically  evaluated  Green's  function  in  an  anisotropic  crystal57  is  presented  in 
Sec.  II.  This  procedure  takes  proper  account  of  the  noncollinearily  of  the  ray  and  propagation 
directions  in  an  anisotropic  crystal.  In  Sec.  Ill  the  Green’s  function  solution  is  used  to  find  an 
expression  for  the  scattered  power  inside  the  crystal.  The  expression  differs  in  several  ways 
from  the  best  previous  treatment  of  Brillouin  scattering  in  anisotropic  media  by  Molulevich. 1 1 
One  difference  is  the  appearance  of  the  Gaussian  curvature  K of  the  surface  w(k)  =oj  in  our 
formula  for  the  scattered  power  inside  the  crystal  Motulevich's  treatment  was  not  based  on  a 
Green’s  function  approach  but  rather  on  Ginzburg’s12  Hamiltonian  approach  using  an  aniso- 
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tropic  Coulomb  gauge. 

In  Secs.  IV  and  V we  present  the  remedy  to  the  second  deficiency  of  previous  theories, 
that  is,  to  relating  the  scattered  power  inside  the  crystal  to  quantities  characterizing  the  detector 
outside  the  crystal.  In  Sec.  IV  the  solid  angle  inside  the  crystal  is  related  to  the  expanded  solid 
angle  outside  the  crystal.  The  formulas  presented  are  completely  general  in  applying  to  any 
orientation  of  the  crystal  axes,  the  scattered  ray,  and  the  surface  normal.  In  Sec  V the  length 
of  the  scattering  volume  along  the  incident  beam  is  related  to  the  corresponding  demagnified 
length  as  seen  outside  the  crystal  by  a field  stop  of  the  detection  optics.  1 he  demagnification 
formula  is  also  completely  general  in  applying  to  any  orientation  of  the  crystal  axes,  the  scat- 
tered ray,  and  the  surface  normal.  The  width  of  the  scattering  volume  (the  unscattered  beam) 
when  laser  sources  are  used  is  typically  much  less  than  the  width  of  the  field  stop  of  the  detec- 
tion optics  and  so  does  not  enter  the  formulas  explicitly  The  solid  angle  expansion  and  source 
volume  demagnification  expressions  tire  then  combined  with  the  scattered  power  formula  of 
Sec  111  to  produce  in  Sec.  VII  a scattered  power  formula  applying  to  measurements  made  out- 
side the  crystal. 

The  scattered  power  formulas  of  Secs  111  and  VII  are  formulated  in  terms  of  an  arbitrary 
mechanism  of  light  scattering.  In  Sec.  VIII  the  power  formula  is  specialized  to  two  important 
mechanisms,  Brillouin  scattering’  ''  (from  acoustic  phonons)  and  Raman  scattering’4  (from 
optic  phonons).  The  final  formulas  are  completely  general  yet  compact  and  convenient. 

II.  INSIDE  GREEN'S  FUNCTION 

The  wave  equation  Tor  the  electric  field  E(r'exp(  —iwt)  m an  anisotropic  dielectric  gen- 
erated by  the  nonlinear  polarization  P N/  ( r ' ) exp  ( — iiui ) at  the  single  frequency  u>  is 

,.2  ,,2  i>  x/ 

V x ( V x E ) — —7-k  ( to  ) E = — -r-  (2.1) 

c2  c2  e„ 

It  has  the  Green’s  function  solution4 
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g*(R)  , f*  -CXP.  !'k*  ‘ ft1  . (2.7) 

47rK 

where 

f*  = -C0S^L  (2.8) 

A*Va* 

and  A'*  is  the  Gaussian  curvature  of  the  w (k ) surface  at  k*,  Equation  (2.7)  has  been  shown7 
to  agree  precisely,  in  the  uniaxial  case,  with  the  Green’s  function  obtained  without  asymptotic 
approximations.  This  Green’s  function  disagrees  by  the  factor./'*  with  the  intuitive  notion  that 
for  a given  direction  of  observation  one  may  use  the  Green’s  function  for  an  isotropic  medium 
with  the  index  of  refraction  appropriate  to  direction  k*. 


III.  SCATTERED  POWER  INSIDE  THE  CRYSTAL 

If  Eqs.  (2.2),  (2.5),  (2.7)  and  (2.8)  arc  used  to  calculate  the  electric  field  and  if  the 
corresponding  magnetic  field,  H = (VxE)/(/w/x(l),  is  also  found,  then  the  Poynting  vector 
for  a given  mode  <t>  is  found  to  be 

s*  = >y»4  »*|/*|2|c*|2t 
32  ir2c  (cosfi* ) 3 

where  //*  is  the  index  of  refraction  appropriate  to  wavevcctor  k'b  and 

e*  = e*  • J P'v,  (r)exp(  — /k*  ■ r)  dr  . 


(3.1) 


(3.2) 


The  ratio  of  the  scattered  power  inside  the  crystal, 

Fins"  -|S*|f2c/nr. 

(where  d fi r is  a solid  angle  of  rays  in  r space)  to  the  incident  power, 


Csc  " A IS")  = ‘cos8H 


may  be  expressed  as 


oscai 
1 ms 

tJ,nc/A 
‘ ins  ' 


= p i/S. 


fil'V/O,; 


(3.3) 

(3.4) 

(3.5) 


Here  F5  is  the  scattering  volume  accepted  by  the  detector.  A,  the  cross-sectional  area  of  the  in- 
cident beam  inside  the  crystal  and  R.  the  scattering  efficiency  (the  scattered  power  per  unit  in- 
cident power,  per  unit  solid  angle,  per  unit  path  length),  is  given  by 


R 


8rr  Vcos8*co.s5h  (**)2A* 

The  nonlinear  phenomena  that  give  rise  to  the  scattering  are  included  in  the  quantity 

J*  ■’ V' 


(3.6) 


J = 


(r)  rdr 


i ' 


. 2 


E" 


(3.7) 


Vs 


One  striking  way  that  the  expression,  Eq.  (3.6),  for  the  scattering  efficiency  differs  from  previ- 
ous expressions  is  by  its  dependence  on  the  Gaussian  curvature  A*  of  the  <t>  branch  of  the 
<o(k)  = (o  surface. 
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IV.  SOLID  ANGLE  EXPANSION 

The  ratio  of  the  solid  angles  subtended  by  the  rays  of  a beam  inside  and  outside  a crystal 
can  be  factored  as 

ds\;n  dnk„ 

‘'“ out  rfHj*  ‘'“out  ' 

The  first  factor  describes  the  ratio  of  the  solid  angles  in  ray  vector  space  and  wavevector  space 
and  is  independent  of  the  existence  of  a crystal  surface.  The  second  factor,  describing  the 
change  in  wavevector  solid  angles,  is  completely  determined  by  Snell’s  law. 

Since  dfl A is  readily  computed  in  terms  of  the  area  of  a patch  dAk  of  the  <u(k)  surface, 
see  Fig.  2,  and  d fl  'n  is  related  to  the  same  area  by  Gauss’  theorema  egregium,  see  Fig.  3,  we  ob- 
tain our  previously  quoted  result7 

dil'Jdn^  = K*(**)2/cos8*  . (4.2) 

Across  the  surface  of  a crystal  (nominally  in  the  3 or  z direction)  Snell’s  law  guarantees 
the  continuity  of  the  transverse  components  of  the  propagation  vector  and,  hence,  of  dk\dk2. 
Since  dk\dk2  is  simply  related  to  the  patch  area  dAk,  see  Fig.  4,  and  hence  to  the  solid  angle 
dil  k,  sec  Fig.  2,  Snell’s  law  leads  to  the  relation 7 

‘'“in  cosS^cosa  ,A 

—7: = — —r . (4.3) 

‘'“oui  (//*) 2 cos/3 

where  /3  is  the  angle  of  arrival  of  the  ray  inside  the  crystal  to  the  surface  normal,  8*  is  the  an- 
gle between  ray  and  wavcvectors,  as  before,  and  a is  the  angle  between  the  departure  ray  out- 
side the  crystal  and  the  surface  normal. 

The  product  of  Eqs.  (4.2)  and  (4.3)  yields  the  desired  solid  angle  expansion 

dfljrn/d(loul  = (ot/c)2 K^cosa/cosp  . (4.4) 

A slight  rearrangement  of  this  equation  suggests  that 

dil  'cos/VA'  (4.5) 

is  an  invariant  for  the  passage  of  a beam  from  one  material  to  another,  a result  we  have  recent- 
ly proved  quite  generally.  Equation  (4.4)  is  a special  case  of  this  invariance  in  which  the 
second  medium  is  a vacuum  with  Avat  = (c/w)2. 


V.  SOURCE  VOLUME  DEMAGNIFICATION 


When  the  laser  beam,  scattered  ray,  surface  normal,  and  departure  ray  outside  the  crystal 
are  all  in  one  plane,  it  is  possible  to  derive  the  length  ls  along  the  laser  beam  in  the  crystal 
from  which  radiation  is  admitted  by  a detector  field  stop  of  length  lD.  As  seen  from  Fig.  5, 
these  lengths  arc  related  by 


/ysindy 


I 


I, 


cos  fi  cos/3  coso 

where  Wy  is  the  scattering  angle. 

When  the  rays  mentioned  above  arc  not  all  coplanar,  Eq.  (5.1)  must  be  replaced  by 


(5.1) 


'.V  = N cosj3 
Id  sinWyCosa 

where  /Vis  the  noncoplanarity  correction, 16 


(5.2) 
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/V  = l/(cos<£  cos </>'  — cos/3  sin<£  sin<£')  . (5.3) 

Here  <t>'  is  the  angle  of  tilt  between  the  plane  defined  by  the  unscatlered  laser  beam  and  the 
normal  to  the  input  surface  and  the  plane  defined  by  the  scattered  ray  and  exit  surface  normal; 
<£  is  the  angle  of  tilt  between  the  latter  plane  and  the  departure  plane  defined  by  the  normal  to 
the  exit  surface  and  the  scattered  ray  after  it  has  left  the  crystal. 

VI.  SCATTERED  POWER  OUTSIDE  THE  CRYSTAL 

To  convert  Eq.  (3.5)  to  a power  scattering  formula  outside  the  crystal,  we  use 

oscai  _ oscal /TCxii  (6.1) 

Ills  1 OUl  ' ' 

pmc  = P int  7«nl  (6.2) 

ms  out 

where  T001  and  rcxil  are  entrance  and  exit  transmission  factors.  If  we  note  that  Vs/ A = /j, 
Eqs.  (4.4)  and  (5.2)  can  be  combined  to  yield 

w 2 Nl,}tl ill>K+  .. 

ls<M{n  = - ~7. (6.3) 

’ 11  c s)r>Os 

where  we  have  written  d ill}  for  di  1,,,..  to  remind  us  that  it  is  the  detector  solid  angle.  A de- 
tailed  derivation  of  this  result  has  been  given  elsewhere.  • 1 • 1 - ■ 1 

If  Eqs.  (6.1),  (6.2)  and  (6.3)  are  combined  with  Eq.  (3.5),  we  obtain 
psca,  Rr'"T»%d(lDNK*  (J2 


If  Eq.  (3.6)  is  used  for  the  scattering  efficiency,  we  obtain 

L 4 MpM”  (6  5) 

P™  \ C sintf v 87r2;t<<>/twcos5<(,cos8H 

Equation  (6.5)  incorporates  all  the  geometric  optics  of  the  crystal  surface.  The  mechanism  of 
the  scattering  is  contained  in  7,  whose  evaluation  will  be  discussed  in  the  next  section. 

VII.  APPLICATION  TO  LIGHT  SCATTERING 

To  apply  our  final  scattering  equation,  Eq.  (6.5),  to  Brillouin,  Raman,  or  some  other 
scattering  mechanism,  it  is  necessary  to  evaluate  J,  of  Eq.  (3.7),  which  involves  the  nonlinear 
optical  properties  of  the  scattering  medium.  Since  the  volume  (/V  is  large  compared  to  all 
relevant  wavelengths,  it  is  permissible  to  take  the  limit  as  Vs  approaches  infinity.  The 
Weiner-Khinchin  theorem17  applied  to  spatial  variables  rather  than  the  lime  then  permits  Eq. 
(3.7)  to  be  rewritten  as  an  autocorrelation, 

ef  ea  J’  exp  ( - r k * ■ r ) < P,NL  (0  )*P?L(r)  > 

J=  2«f?  |E"(r)|2  ’ 

where  the  limits  arc  infinite. 

With  the  understanding 

E"(r.r)  = y [ Ew(r ) + E" (r ) * em"  ] . (7.2) 

PNL(t .,)  = | [p  NL  (r ) e ■'"*'  + P vt  (r ) * e**]  . (7.3) 

where  o> L is  the  input  or  laser  frequency  and  wp  is  the  Brillouin  (or  Raman)  scattered  frequen- 
cy, we  can  write 


'U+  **' 
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C'  <r./>  = *„  (r./)  Ztf  (r./>  (7.4) 

where  the  "susceptibility"  *„/,  is  space  and  time  dependent  because  it  is  induced  by  the  pho- 
non field  in  Brillouin  or  Raman  scattering.  Equation  (7.1)  can  then  be  simplified  to 

J =4-  ( exp  [— /(k*-kw)  rl  < N**i0.0)  N*"l r.O)  > dr  (7.5) 

2 J 

where 

N*"  (r.O  = ef  Xah  (r.O  ^ (7.6) 

and  we  have  used  the  plane-wave  character  exp  ( / k H • r ) of  the  unscattered  wave  Ew(r). 

If  one  is  concerned  with  line  shape,  Eq.  (7.5)  can  be  decomposed  by 

oo 

J — f J( oO  do>  /2tt  (7.7) 

Jo 

where 

y(w)  = yJJ  exp  (/(w  — atL  )r]  exp  (-/(k* -kH ) • r]  < (0.0 ) (r.O  > dr  di,  (7.8) 

and  only  the  positive  frequency  components  of  exp  (—iotL  i)  N (r.O  are  included  so  that 
J(ot)  vanishes  for  ot  < 0. 


Brillouin  Scattering 

The  first  application  of  our  modified  Green’s  function,  Eq.  (2.7),  with  surface  corrections 
was  made  in  connection  with  a detailed  study  of  Brillouin  scattering  in  calcite.  Because  the 
formulas  in  this  paper  give  a clear  factorization  of  our  Brillouin  scattering  formula1  into  intrin- 
sic and  geometric  components,  we  indicate  here  the  evaluation  of  J based  on  Eq.  (7.5).  To  be 
consistent  with  our  definition  1 of  the  photoelastic  susceptibility  as  a relation  between  the  posi- 
tive frequency  components  of  P*L,  EH  and  ucd  (the  displacement  gradient),  we  must  write 

Xah  ~ 7 Xahcd  uc.J  ' (7.9) 


The  averages,  < u(i  (0.0)  u t (r ,/) 
terms  of  normal  coordinates,15 


> , can  be  evaluated  by  using  the  expansion  of  uc  in 


u (r.z ) = £ 

q 

where  q = k'1  is  the  acoustic  phonon  propagation  vector,  p is  the  crystal  density  and  fl  is  its 
volume  so  that  p 11  = MN  = mass  per  unit  cell  x number  of  cells.  The  unit  displacement  vec- 
tor b is  characteristic  of  the  type  of  mode  (e  g.  transverse),  and  the  sum  over  q also  implies  a 
sum  over  types  of  modes.  The  mode  amplitudes  in  the  quantum  mechanical  case  obey 

< af(q.i)  a(q',i)  > =/75(q,q'), 

< o(q,/)  af(q',i)  > =(«+!)  8(q,q') 

where  ii  \ s the  actual  phonon  excitation  number  that  reduces  to 

- j 

exp  [1}ot(q)/kT\  —1 

in  the  thermal  equilibrium  case.  Averages  of  the  type,  Eq.  (7.12),  contribute  to  Stokes  scatter- 
ing, whereas  those  of  type,  Eq.  (7.1 1),  contribute  to  anti-Stokes  scattering.  Because  we  usually 
have  1\a>«kT, \ both  Stokes  and  anti-Stokes  scattering  have  an  intensity  proportional 
ii  = kT/Tiot.  If  we  combine  Eqs.  (7.5)  and  (7.9M7.I3),  we  obtain 


(7.11) 

(7.12) 

(7.13) 


ft 

2 p n at  (q  ) 


b e'q  r « (q.r ) + e iil  Tat(q,i)  . (7.10) 
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j = — - G , C=|  ef  ef  Xjjki  bk  ai  1 2 (7,14) 

where  a,  = q ,/  |q|  = kj4/  | k | is  a unit  vector  in  the  direction  of  phonon  propagation  and 

vA=ai(kA)/kA  (7.15) 

is  the  sound  velocity. 

Equation  (6.5)  for  the  Stokes  (or  anti-Stokes)  power  scattering  by  a single  type  of  acoustic 
mode  can  be  rewritten  as 


bsc  a i 
1 out 

pine 
1 OUl 


4 niL)  m 0 


rmrxil  k tg 


sin«s-  8rr2  /;H  cosS4  cos8wpvj 


(7.16) 


where  O of  Eq.  (7.14)  can  also  be  expressed  in  terms  of  the  Pockels  tensor,  pijk/,  by  means  of 

C = (1/4)  (n’t>  n")4  (cos8  * cos8  H ) 2 F , (7.17) 

F=  | df  d"  plJk,  bk .a,  \2  . (7.18) 

Here  dH  and  are  unit  electric  displacement  vectors  of  the  input  and  scattered  beams,  and  o> 
is  the  Brillouin  scattered  frequency  id14  . Equation  (7.16)  was  used  in  the  analysis  of  our  exper- 
imental results1,19  except  for  the  noncoplanarily  factor.  A/,  which  was  not  needed  then  since  all 
experiments  were  conducted  in  a symmetry  plane. 


Raman  Scattering 

A detailed  analysis  of  Raman  scattering  by  polaritons  using  the  fluctuation-dissipation  ap- 
proach of  Barker  and  Loudon6  has  been  given. 14  To  evaluate  Eq.  (7.8)  for  Jiui),  needed  to  ob- 
tain the  line  shape  in  Raman  scattering,  we  express  the  nonlinear  susceptibility  defined  in  Eq. 
(7.4)  in  the  usual  way. 


X„  (r./)  = 2£ 


AS 


vM  (r.t)  + 2 Bjjk  Ek  (r,/) 


(7.19) 


as  an  ionic  contribution  associated  with  the  displacement  >vM  of  mode  p plus  an  electronic  con- 
tribution proportional  to  the  electric  field  £(r,r).  The  factors  of  2 assure  conformity  with 
definitions  in  our  previous  work. 14  Because  field  and  displacement  are  correlated  in  polariton 
motion,  it  was  simpler  to  evaluate  the  correlations  with  the  help  of  the  fluctuation-dissipation 
theorem. 14,6  If  Ts.  represents  the  response  of  < B>  to  a unit  force  at  frequency  w applied  to 
A,  Eq.  (7.8)  yields  ■20 


7(w  + (i>L)  = 41)  n(ai)  cf  e”  e$  e"  lm[  Jijah  ] 


(7.20) 


where 


Ji/ab 


= 1 ^V'Aab  V + TcC 


+ I Up  > ’ Babc  TCII,  + I B’k  A'ah  Tvk,  (7.21) 


I KApV  «"(«)  + B’k 


V 


LASb  < + Bahc 


+ E up  ) * Aph  B»  (ai)  + £ {Apr  B>L  Uo ) 


(7.22) 


The  second  form,  Eq.  (7.22),  has  used  the  equations  of  motion  to  express  all  response  func- 
tions in  terms  of  the  field-field  response,  i.c.  the  response  of  the  electric  field  to  a unit  applied 
external  polarization,  as 
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rkt  = Ejn"  = («_l  >,a/« 

= I 


s..  S, 


*-1.2 


t cos2  8*  [(t*/o»)2  - (//^  ) 2 ] 


(7.23) 


e„  s • 


K (to  ) 


• s 


The  first  two  terms  in  Eg.  (7.23)  represent  the  transverse  modes  while  the  third  term 
represents  t lie  longitudinal  mode. 

The  mode  w**  will  generally  obey  an  equation  of  the  form20 

m"  [ (w"  )2  - o>2  - iio  (to ) | tv**  = q?  Ec  (7.24) 


where  nt>L  is  the  elfeclive  mass,  the  elfective  (angular)  frequency,  y^  (at)  the  effective  fre- 
quency dependent  damping  constant  and  qf  is  the  c component  of  the  effective  charge.  The 
response  of  to  a unit  applied  force,  needed  in  Eq.  (7.22),  is  given  by 


ir  (to)  = 


nr 


IU"); 


— (O  2 — 


iw  y M (a>  ) ] 


(7.25) 


We  may  combine  Eqs.  (7.20)  — (7.25)  to  obtain  J(to  + to2-).  If  V(ta  + tnL)  is  inserted 
in  place  of  ,/  in  Eq.  (6.5),  the  result,  in  view  of  Eq.  (7.7),  is  the  ratio  of  (he  Raman  scattered 
power  per  unit  frequency  interval  to  the  incident  power,  both  computed  outside  the  crystal. 
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FIGURE  CAPTIONS 

1.  A typical  experimental  setup  for  a Raman  scattering  experiment  which  displays  the 
expansion  of  the  solid  angle  on  emerging  from  the  crystal.  For  a narrow  laser 
beam,  the  source  volume  I s — Als  from  which  scattered  light  is  accepted  is  limited 
by  the  input  beam  area  A and  a length  ls  determined  by  the  field  stop. 

2.  The  solid  angle  dfl*  associated  with  a patch  of  area  dAk  — dudv  on  the  surface 
to(k)  =o»  of  free-wave  k vectors  is  given  by  difk=dAk  cos8/A2,  where  dA^cosS  is 
the  component  of  the  area  dAk  normal  to  k,  since  dAk  is  normal  to  t or  Vto  (k). 

3.  Gauss's  bump  theorem  relates  the  area  dA  = dudv  of  an  element  of  surface  to  the 
area  d if.  ’ on  the  surface  of  a unit  sphere  subtended  by  the  unit  normals  to  dudv 
that  have  been  shifted  in  a parallel  manner  until  their  origins  coincide  at  0,  the 
sphere  center.  The  solid  angle  dif  = dh  d<t>  where  dl>  — du/p„  is  the  angle  AUD 
and  d<l>  = dv/py  is  the  angle  AVB  where  pu  and  p are  the  principal  radii  of  curva- 
ture. Thus  di\r  — K dAk  where  K = ( pu  ) “'  is  the  Gaussian  curvature  and 
dAk  = dudv  is  the  area  of  the  patch  on  the  m (k  ) surface. 

4.  The  patch  dAk  = dudv  on  the  surface  cutk'  = w of  Fig.  3 is  plotted  in  a Aj  ,k2 .Aj 
coordinate  system  in  which  the  three-direction  is  along  the  direction  n of  a unit 
normal  to  the  surface  of  the  crystal.  The  projection  of  dAk  onto  the  1-2  plane 
yields  dA,  dk2  — dAk  cos/3  because  £ is  the  angle  between  the  normal  t to  '.he 
patch  (the  ray  direction)  and  the  normal  n to  the  surface. 

5.  Demagnification  corrections  when  the  arrival  ray,  departure  ray,  and  surface  normal 
are  all  in  one  plane.  The  detector  field  stop  (sec  Fig.  1)  is  represented  in  image 
space  by  the  knife  edges,  which  accept  a dimension  lt)  perpendicular  to  the  beam. 
The  portion  ls  of  the  laser  beam  accepted  is  determined  by  the  geometrical  condi- 
tions shown,  independent  of  the  orientation  of  the  virtual  image 
lA  : /ssin«s/cos/3  = l/cosfi  = Ip/cosa. 


This  manuscript  was  produced  on  a computer  driven  phototvpesetter  using  the 
UNIX  operating  system  developed  at  Bcli  Laboratories. 
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